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THROUGH NEUTRINO OBSERVATIONS 
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Abstract. The origin of ultra-high energy cosmic rays (UHECR) has 
been widely regarded as one of the major questions in the frontiers 
of particle astrophysics. Gamma ray bursts (GRB), the most violent 
explosions in the universe second only to the Big Bang, have been a 
popular candidate site for UHECR productions. The recent IceCube 
report on the non-observation of GRB induced neutrinos therefore at- 
tracts wide attention. This dilemma requires a resolution: either the 
assumption of GRB as UHECR accelerator is to be abandoned or the 
expected GRB induced neutrino yield was wrong. It has been pointed 
out that IceCube has overestimated the neutrino flux at GRB site by a 
factor of ~ 5. In this paper we point out that, in addition to the issue 
of neutrino production at source, the neutrino oscillation and the pos- 
sible neutrino decay during their flight from GRB to Earth should fur- 
ther reduce the detectability of IceCube, which is most sensitive to the 
muon-neutrino flavor as far as point-source identification is concerned. 
Specifically, neutrino oscillation will reduce the muon-neutrino flavor 
ratio from 2/3 per neutrino at GRB source to 1/3 on Earth, while neu- 
trino decay, if exists and under the assumption of normal hierarchy of 
mass eigenstates, would result in a further reduction of muon-neutrino 
ratio to 1/8. With these in mind, we note that there have been efforts 
in recent years in pursuing other type of neutrino telescopes based on 
Askaryan effect, which can in principle observe and distinguish all three 
flavors with comparable sensitivities. Such new approach may therefore 
be complementary to IceCube in shedding more lights on this cosmic 
accelerator question. 
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1 Introduction - Neutrino as Cosmic Messenger 

Ultra-high energy cosmic rays (UHECR) with energies beyond 10 19 eV have been 
observed (HiRes 120091 Auger l20l0l ANITA |20T0|) . Their origin, however, has 
not been clear. This challenge has been considered as one of the eleven science 
questions for the new century by the well publicized white-paper prepared for the 
US National Research Council by the Turner Committee in 2003 (Turner 2003). 
Gamma ray bursts (GRB), the most violent explosions in the universe second 
only to the Big Bang, has long been considered as a promising candidate site 
for the "cosmic accelerator" where UHECRs, primarily protons, are generated. 
Since proton trajectory can be bent by inter- and intra- galactic magnetic fields 
on its way to Earth, its incoming angle cannot reveal its origin. It happens that 
a necessary by-product of UHECR on or near the cosmic accelerator, such as 
GRB, site are the neutrinos at comparable energy (UHECN). Being charge- neutral 
and weakly interacting, neutrinos so produced can propagate straight-forwardly 
to Earth. The detection of such neutrinos would therefore provide a useful means 
to address the cosmic accelerator puzzle. 

There are multiple channels where such conversion can happen. We assume 
that the UHECRs are primarily protons. First, such protons can interact with the 
intense GRB background within the GRB fireball, possibly comoving, and neutri- 
nos are generated through photo-pion production. Such interaction is particularly 
pronounced in the A-resonance process: 

p + j^A^n + n + ^n + + v^} ->• n + {[e + + v e + v^] + v^}. (1.1) 

The resulting UHECN spectrum should peak at several hundred TeV (Waxman 
& Bahcall [1997|) . Such UHECNs are produced strictly on the GRB site. Second, 
the UHECR protons may instead collide with the cosmic microwave background 
(CMB) on their flight to Earth, and turn into UHECNs under the same process 
in Eq.(l). This is the so-called GZK process (Greisen 119661 Zatsepin & Kuzmin 
119661 Berezinsky & Zatsepin I1969[) . For this channel to happen, it would require 
a much higher UHECR proton energy since CMB photons are much softer than 
that of GRB (but may collide head-on). The associated neutrino spectrum has 
a flat top ranging from 10 15 eV to 10 20 eV (see Fig. 1). Though such UHECNs 
are produced offsite from the cosmic accelerator, it generally occur very near the 
source. For example, at redshift z ~ 1 the mean-free-path for a UHECR proton 
above the GZK threshold energy is < 6Mpc. A neutrino so produced would arrive 
at Earth with an incoming angle that is within 5 arc second around the GRB 
where the UHECR proton was accelerated. The point-back ability of UHECN to 
the cosmic accelerator is thus not so compromised even through this channel. 

2 Non-observation of GRB Neutrinos at IceCube 

Giving the importance of the issue as described above, the recent IceCube an- 
nouncement of the non-observation of GRB neutrinos (Abbasi et al. [IceCube 
Coll.] 120111 Ahlers, Gonzalez-Garcia & Halzen 2011} came as a disappointment 
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to many. IceCube calculated the expected prompt neutrino spectrum in the in- 
ternal shock scenario of the fireball model following Guetta et al. {2004), which is 
based on Waxman & Bacall (|1997j) . As far as point-source observations is concern, 
IceCube is most sensitive to muons for two reasons (Karle et al. [IceCube Coll.] 
2003): 1.) Muons allow a very good angular resolution of 0.7° over a wide range 
of energy; 2.) The effective volume for muons exceeds the geometric volume of the 
detector by factors of 10-50, depending on energy. Following IceCube and for the 
rest of this paper, we shall neglect the contributions from the other two flavors in 
IceCube's sensitivity to address the GRB connection with cosmic accelerator. 

IceCube's starting point is to express the GRB neutrino flux in terms of its 
UHECR proton flux: 

F lc 1 

J-p 8 

where f^.b = f-n{E = Eb) is the fraction of proton energy carried by the pion 
as a result of proton-GRB photon collision, with a spectral-break energy at 
A-resonance: 

E b = 1.3 x 10 16 r| 5 e^ oV eV, (2.2) 

where T = 10 2 5 r2.5 is the Lorentz factor of the bulk flow and eb = lefc.McV- On 
the other hand, the proton flux can be normalized by the gamma-ray flux through 



J-p — J-ry, (2-3) 

Je 



where f e is the ratio of accelerated proton to electron energies. By 2011 with half 
of its detector completed, IceCube has reached the sensitivity that is comparable 
to the expected neutrino flux from GRBs. Based on 117 GRB events, IceCube 
searched for neutrino emission from these sources over a wide range of energies and 
emission times, but produced no evidence for such, excluding prevailing models at 
90% CL. 



3 Over-estimation of Neutrino Flux at Source 



It has been pointed out, however, that the IceCube Collaboration might have 
over-estimated the neutrino flux at GRB source by as much as a factor 5 (Zhuo Li 
2012). If so, then the IceCube results have not yet ruled out the GRB candidacy 
for the cosmic accelerator. Here we review Li's argument and correction to the 
GRB neutrino flux. His main points are as follows. 

1. For a flat proton distribution with the index given above and the Band- 
function parameters, a = — 1 and /? = —2.2, so chosen by IceCube, the assumption 
of f n (E) — f Wt b is valid only for protons with energy E > E\>. For E < Eb, Jv is no 
longer a constant but instead proportional to the energy: jv oc E. The physical 
reason for this is simple: the higher energy target photons are harder to find to 
match with the lower energy protons. Thus by applying constant f n to all proton 
energies, Eq.(l) has overestimated the initial neutrino flux. 
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2. Eq.(3.1) has also ignored the suppression of neutrino production at high 
energies due to the radiative cooling of secondary pions/muons. The synchrotron 
cooling timescale is shorter than that for the secondary decay, and therefore the 
neutrino production will be suppressed, when the pion/muon energy is above the 
cooling energy E c . 

Thus the neutrino production is mainly contributed by the primary protons 
with an energy window Ef, < E < E c , which is only a fraction of the total number 
of GRB accelerated UHE protons. Let the maximum and minimum accelerated 
proton energies be E max and E min , respectively, with ratio E max / E min ~ 10 9 , and 
let E c /Eb ~ 10 2 . One then finds, for a proton energy spectrum E 2 dn p /dE oc E 2 ~ p , 
with p « 2, 



Thus the correction to the IceCube assumption of the neutrino flux is the reduction 
of Fl c by roughly a factor 5. With this correction, Li found that the IceCube data 
is still consistent with the assumption of GRB as the origin for UHECR. 

It should be reminded that there exist additional corrections to the neutrino 
flux due to various uncertainties in the fundamental physics parameters such as 
f„ and f e , as well as in the GRB and UHECR acceleration models. 

4 Evolution of Neutrino Flavors in Flight 

Having discussed the conversion factor from UHE protons to neutrinos at GRB 
site, which is GRB model and cosmic accelerator model dependent, in this sec- 
tion we turn our attention to the impact on the GRB-neutrino connection due to 
neutrino properties in-flight from GRB to Earth. While the notion of neutrino 
oscillation is by now well-known and has been incorporated into data analysis, 
it appears that the implication of possible neutrino decay has thus far not been 
taken into active consideration. 

4.1 Neutrino Oscillation 

As we know, neutrinos oscillate among their three different flavors (Fukuda et al. 
[Super-Kamiokande Collaboration] 1998). From Eq.(2.1) we see that every proton- 
photon interaction will produce one electron neutrino and two muon neutrinos. So 
the ratio between different neutrino flavors at source, per out-coming neutrino, is 



Being produced at cosmic distance, which is much much larger than the neutrino 
oscillation length, the GRB induced ultra-high energy neutrinos flavors will reach 
an equilibrium when arriving at Earth with a uniform distribution (see, for exam- 
ple, Wang, Chen, Huang et al. [2013 ): 




(3.1) 



At Source : 



/f : J* : / r s = 1/3 : 2/3 : 0. 



(4.1) 



In Flight : 



./. : : fr = V3 : 1/3 : 1/3. 



(4.2) 
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As mentioned earlier, although IceCube is able to detect all neutrino flavors, in 
its search for point courses such as neutrino emission from GRB, muon- neutrino 
has the highest sensitivity and best angular resolution. This, however, necessarily 
reduces the effective neutrino flux to f/3 of the total flux emitted from GRB. 

4.2 Neutrino Decay 

Another major discovery about neutrinos in the last two decades is that neutrinos 
have mass. This, together with additional assumption of the decay process based 
on notions beyond the standard model of particle physics, leads to the prediction 
of neutrino decay, from the higher mass eigenstates to the mass ground state. 
One can envision at least two possible mass hierarchies: the normal hierarchy, 
i.e., to 3 3> m 2 > mi, and the inverted hierarchy, i.e., m 2 > mi 3> m 3 . It 
has been shown that under normal hierarchy the eventual flavor ratio on Earth, 
again normalized to number of neutrinos produced at GRB site, is (Beacom, Bell, 
Hooper, et al. |2"UU51 Maltoni & Winter H]® 

Normal Hierarchy : ff : : /f = 2/3 : 1/8 : 5/24, (4.3) 

or 

Inverted Hierarchy : /f : f® : jf = : 2/5 : 3/5. (4.4) 

So for the IceCube search for GRB emission of neutrinos, which is most sensitive 
to its detector sensitivity would be further reduced from 1/3 (based on pure 
oscillation) to 1/8 of the total neutrino flux for the case of normal hierarchy, 
but would gain slightly from 1/3 for the purely oscillating scenario to 2/5 if the 
neutrino mass eigenstates follow the inverted hierarchy. 

In summary, at the GRB site every UHECR proton that collides with the 
background GRB or CMB photons would generate 3 neutrinos, among them 2 
are muon neutrinos, i.e., the contributes 2/3 of the total population. By the 
time these neutrinos arrive on Earth, the contribution reduces to 1/3, due to 
neutrino oscillations. When the neutrino decay is considered, in particular with 
the assumption of normal hierarchy then its contribution will become 1 /8 when 
arriving on Earth, a factor 3/8 reduction from that under pure oscillation. On the 
contrary, if the neutrino mass eigenstates follow the inverted hierarchy, then the 
ffx population rises from 1/3 to 2/5. We should caution, however, that the notion 
of neutrino decay is not a direct consequence of the standard model and is not yet 
experimentally verified. 

5 Detecting GZK Neutrinos with Askaryan Effect 

While the estimate of the GRB induced UHE proton and neutrino fluxes at the 
GRB site are GRB model and acceleration mechanism dependent, we see that 
the evolution of the neutrino flavor ratio during their flight to Earth depends 
on neutrino's fundamental properties. Among them neutrino oscillation is well- 
established, while neutrino decay still awaits experimental evidence. Based on the 
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discussion in the previous section, however, it should be fair to conclude that, 
as far as the issue of UHECR cosmic accelerator is concerned, it would be more 
advantageous if a neutrino observatory can be sensitive to all flavors. In this 
regard, an idea proposed by Askaryan in the 1960s ( Askarvan 119621 [T565)), can 
in principle detect all neutrino flavors with comparable sensitivities due to the 
vastness of the typical detection volume involved in such approach. 

As was proposed by Askaryan, high energy cosmic neutrino can be detected by 
observing the radiowave band of the Cherenkov radiation emitted by the neutrino- 
induced shower in a large solid target. The neutrino induced shower, though 
inherently charge neutral, will develop a charge disparity by the time when the 
shower reaches its maximum, with ~ 20% more electrons than positrons due to the 
shorter stopping distance of the positrons propagating in ordinary matter. Since 
all particles travels near the speed of light, the shower remains compact and thus 
the emitted Cherenkov radiation is a sharp impulse and therefore is wide-band 
in the frequency domain, where the radiaowave portion is enhanced due to the 
coherence of emission. This effect has been validated in a series of experiments in 
sand, salt and ice performed at SLAC (Gorham 2000). 

Several ongoing and proposed experiments, e.g., GLUE (Gorham et al. I2004p . 
RICE (Kravchenko et al. 12006]) . ANITA ([2009"1> . LUNASKA (James et al. l20T0j) . 
ARIANNA (Barwick [20011 Gerhardt et al. l20T0j) . and ARA (Chen & Hoffman 
120091 Allison et al. [ARA Coll.] I201ip . are based on the Askaryan effect. In 
particular, the balloon-borne ANITA has so far completed two successful missions 
in Antarctica with exciting results. With an average altitude of 30 km, it tends 
to detect neutrinos at the higher energy end of the GZK spectrum. So far one 
candidate neutrino event has been identified from ANITA-2 (Gorham et al. I2010|) . 
An upgraded ANITA-3 has been scheduled to launch in December 2013 with the 
expectation of detecting multiple neutrinos. Two new ground-based projects, AR- 
IANNA and ARA, are both at the proof-of-principle stage. For the case of ARA, 
it is envisioned to cover 200km 2 area at South Pole with 37 antenna stations. To 
date (Jan. 2013) 3 stations have been successfully deployed about 200m under ice 
at Pole. Figure 1 shows the projected ARA37 detector sensitivity in comparison 
with other projects. 

Focusing on the issue of GRB as a candidate site for UHECR production and 
its connection to UHE cosmic neutrino detection, what is important for a neutrino 
telescope is the angular resolution and the total flux of incoming neutrinos; the 
composition of neutrino flavors within the total flux does not a priori matter. 
That said, it would be advantageous if a neutrino telescope for this purpose has 
indeed the capability of distinguishing neutrino flavors. IceCube can distinguish 
all three flavors based on different characteristics of their tracks. What about 
Askaryan effect-based observatories such as ARA? Recent investigations indicate 
that flavor identification in ARA37 is in principle feasible (Chen, C.-C, Chen, P., 
Hu et al. I2HT31 Wang, Chen, Huang et al. |2UT3|) . On the other hand, due to 
the vastness of the detector's geometric volume and sparseness of antenna stations 
(2km separation), as a trade-off for its higher sensitivity the angular resolution of 
Askaryan effect-based observatories tend not to be as good as that for IceCube. 
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Fig. 1. ARA37 sensitivity in comparison with other contending projects plotted against 
different GZK neutrino models. 

For example the angular resolution for ARA37 is ~ 6° , which is about an order of 
magnitude worse than that of IceCube muon-neutrino detection. 

6 Conclusion 

We have shown that in addition to the over-estimation by IceCube of the GRB 
included neutrino flux at-source by a factor ~5 as pointed out by Zhou Li, there 
are additional impacts to the neutrino flux in-flight as they traverse the cosmos. 
While the effect due to neutrino oscillations is well-known and has been incorpo- 
rated into analysis, that due to a possible neutrino decay has received less atten- 
tion. If neutrinos do decay and their mass eigenstates follow the normal hierarchy, 
then the net effect is that the flux for muon-neutrino would further reduce from 
1/3 under the pure oscillation scenario to 1/8, which would be quite a sizable 
impact on the IceCube detector sensitivity, assuming that IceCube solely relies on 
muon-neutrinos in its investigation of the cosmic accelerator question. There are, 
however, new types of neutrino observatories based on Askaryan effect that are less 
sensitive to neutrino flavors and would therefore provide almost one order of mag- 
nitude improvement in neutrino detection sensitivity. As a trade-off, their angular 
resolution, on the other hand, will not be as good as that with muon-neutrino at 
IceCube. If such detector can in addition distinguish different flavors, then the 
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measured flavor ratio would provide a crucial information on the nature of UHECR 
production, or the inner- workings of cosmic accelerator, at source such as GRB. In 
conclusion the new type of neutrino observatories based on Askaryan effect maybe 
complementary to IceCube in shedding lights on this very acute question about 
the origin of cosmic accelerators. 
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